Oscillatory behavior of the in-medium interparticle potential in hot gauge system 
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We investigate the in-medium interparticle potential of hot gauge system with bound states by 
employing the QED and scalar QED coupling. At finite temperature an oscillatory behavior of the 
potential has been found as well as its variation in terms of different free parameters. We expect the 
competition among the parameters will lead to an appropriate interparticle potential which could 
be extended to discuss the fluid properties of QGP with scalar bound states. 
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I. INTRODUCTION 



When concerning a hot/dense plasma, Debye screening picture is widely accepted where the in- medium screening 
potential is exponentially damping with distance. The damping rate, namely the so-called Debye mass, is the reciprocal 
of the screening length which scales the effective interactive radius of in-medium particles. At high temperature, the 
mean free path of particle is much larger than the screening length, therefore such kind of system can be treated 
y—i • as weakly coupled gas composed of quasi-particles with effective radius. Nevertheless this short-distance dominant 
£> | Debye screening is not the whole story. Due to the dissimilar analytic structures of the effective boson propagator, 
■ the screening potential may behave distinctly at different situations and approximations, in which large distance 
oscillations become equally significant. 

To make it more pellucid, we follow notations in Ref. 1] and defined the in-medium interparticle potential V(r) in 
' linear response framework as[2| 

o \ v{r) = 9iQi lm r dq _J^ (1) 

P-i 

where r is the distance between two particles, Q\ and Q2 are their charges, Q — (u), q) is the four-momentum transfer 
Qh. with q = |q|, and F is the longitudinal boson self-energy defined by F — (Q 2 /q 2 )Hoo- Obviously, the behavior of 
V(r) is totally determined by the analytic structures of the denominator of the integrand on the complex plane. The 
denominator and the so-called dielectric function are almost the same except for a difference of Q 2 . So we can simply 
say that the form of dielectric function determines the interactive potential. 

Take the Friedel oscillations and Debye screening as examples. In the extremely high density but zero temperature 
environment, the explicit expressions of F for nuclear matter^, 0] and gauge plasmaQ have been obtained by 
ignoring the fermion masses. Due to the existence of sharp fermi surface, the dielectric function contains branch 
cuts, which dominate the contribution of contour integration at large distance and lead the interparticle potential 
oscillate considerably instead of monotonic damping. This phenomenon induced by the sharp fermi surface is named 
as Friedel oscillations [H0,0,[5[. However at non-zero temperature the Friedel oscillations will be strongly suppressed 
because of the smeared Fermi surface at finite temperature. Fortunately, although the branch cut contribution is 
suppressed at finite temperature, the pole contribution stands out, resulting in the so-called Yukawa oscillations P, 0] 
which dominate in a wide range of distance. On the contrary, at high temperature but zero chemical potential, there 
is only pole but no branch cut contribution. At hard thermal loop(HTL) approximation, the poles of the integrand 
are on the imaginary axis@ which grants Debye screening. While in a not very high temperature environment, the 
HTL approximation might not be eligible and thus one may try to give up this approximation and employ a complete 
one-loop self-energy. However, the analytic structures of the self-energy are much more complicated in one-loop order. 
Actually, the position of the pole contains both imaginary and real parts on the complex plane which can only be 
identified numerically. The emergence of the real part of the pole means the interparticle potential is not monotonically 
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FIG. 1: One-loop boson self-energy. The solid line is for the fermion and the double solid lines are for the scalar bound state. 

damping but oscillating instead. For example, if there is a pole in the integrand of Eq.([T]), whose imaginary and real 
parts are defined as qi and q r respectively, the potential can be obtained by performing the contour integral in the 
upper plane and recast as 

V ( r ) = r~2~T"Tir\ ~ i a cos (<lrr) + b sm(q r r)} , (2) 

7r(cr + tr) r 

where a and b are defined as the real and imaginary parts of the residue, 

(q 2 + F(q))' 

q 

with the prime denoting for d/dq. Considering the real part of F is an even function, we can only concern the residue 
in the first quadrant and double it in the final result. From mathematical points of view, when the pole is on the 
imaginary axis, the oscillating factor in the blanket of Eq.@ vanishes as a result of q r = 0, which leads to Debye 
screening potential. In other words, if the pole deviates from the imaginary axis, there exit oscillations. 

Besides the mathematical description of the oscillations, we have an another potential motivation on the discussion 
of the oscillatory potential. As we all know the quark-gluon plasma (QGP) produced at rclativistic heavy ion collider 
(RHIC) is almost likely to be perfect fluid at temperatures in the order of several times of critical temperature T ~ (1 — 
3)T C . The fascinating low viscous mechanism has attracted much attention but still challenges both experimentalists 
and theoreticians. One possible way to understand it is to calculate the shear viscous coefficient. In the weakly- 
coupled regime many publications such as Refs. 0i,i,[iElliIiiIiElilIiEIlEISIlElISS Ej HI 0, H El , no 
matter from kinetics or Kubo formula, are overestimating the shear viscosity extracted from experiments or lattice 
results [H,|23, Hl| . The failure prompts to examine the picture of weak coupling and short-distance screening potential, 
and consider the strong coupling and/or long-distance correlation. Apart from the attempts in weak coupling, some 
scientists proposed strongly-coupled plasma and the liquid state theory[29|, [3(| to explain the perfect fluid QGP. One 
should keep in mind that the necessary condition to form a liquid state is the interparticle potential is strong enough 
in coupling as well as long enough in distance [3l|, Obviously Debye potential is not a good candidate because of 
its short distance domination character. Therefore it makes sense to study the oscillatory potential with respect to 
all kinds of dynamical and thermodynamical parameters so as to make sure if it is possible to produce a liquid state 
in such a hot environment. 

Shuryak and Zahed[33j raised the idea of coloered loosely binary bound states to explain the small viscosity at the 
near critical temperatures. In their picture, quark and anti-quark are bounded together, leaving a weak Coulomb 
potential as a remnant for interacting with the colored particles, such as quarks, gluons or other bound states in QGP. 
They argued the possible existing low- lying resonances led to large "unitarity limit" scattering cross section and 
naturally were expected to play an important role in the transport coefficients. In this paper, we adopt this picture 
of relativistic plasma with bound states, attempting to study their effect on the in-medium interparticle potential 
by using QED and scalar QED (sQED) as toy models. Since we do not intend to describe the real electromagnetic 
system but just quote their interaction forms to see the qualitative properties of the in-medium interparticle potential, 
so the temperature, the coupling strength and the masses of fermion and scalar bound state, are free parameters in 
the following calculation. 




ib, 



(3) 



II. QED AND SQED SELF-ENERGY 

In Eq.([T]), one can see all the complication comes from the longitudinal self-energy of gauge boson. In one- loop 
order the photon self-energy is diagrammatically denoted by Fig. 1(a) in pure QED plasma and modified by adding 
FigsQJb) and (c) according to sQED Feynman rules when considering the existence of scalar bound states. The explicit 
expression of F for Fig. 1(a) is composed of vacuum and matter contributions, that are F£ ac and F® natt respectively. 
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The vacuum contribution can be obtained by dimensional regulation and MS scheme. The result can be found in 
many textbooks on Quantum Field Theory [34[, 

where a = is the fine structure constant of QED, and m is the mass of fermion. The matter contribution can be 
obtained either in real time or imaginary time formalism which reads[35| 



Kau(0,<l) = — / dp P 
™ Jo 



Alu a —q, ( q — 2p 
log 



4p q \q + 2p 



»/H. (5) 



where ui is the on-shell energy of fermion constrained to ui = \/ p 2 + m 2 and n,f(u) — (eP^ + is the Fermi-Dirac 
distribution function with j3 = 1/T. 

The existing scalar bound states will " supplement" the polarization of the gauge field and consequently remodel the 
interaction among particles. We involve this effect by employing scalar QED. The polarization and tad pole diagrams 
presented in Figs. 1(b) and (c) will be added to (a). In the real time formalism F can be written as ||36lT37j 

F 6 (0, q) = -2nra s J ^^(2po - k ) 2 [A F (P - K) A R (P) + A A (P - K) A F (P)], (6) 



F c (0, q) = 4i7ra s J ^[A R (P) + A A (P) + A F (P)], (7) 

where a s is the effective strength for the coupling with scalar bound states, capital letters denote for four-momentum, 
Ap, Aji and are the propagators in Keldysh representation[38|, [39( which read, 

(P) 



P 2 — m 2 ± isgn(po)e ' 
A F (P) = -2iri[l + n b (\ Po \)5(P 2 - m% (8) 

where m s is the scalar mass, e is the infinitesimal, sgn denotes for the sign function and rib(\po\) — (e^ ' — 1) _1 is 
the Bose-Einstein distribution function. 

Separating the vacuum contribution from Eqs.([B]) and J7|), one can obtain 



r-l 

47r y J y ' 6 \m 2 + x(l — x)q 2 J 



*l + c c (0, g) - -^q 2 I dx(l - x) 2 log ( ; — ) , (9) 



ittM—— / dpi- 



— log 

p q 



q-2p 
q + 2p 



(10) 



where the vacuum contribution is manipulated in the same renormalization scheme as Eq. @ , where lu s — \Jp 2 -\- ' 



III. NUMERICAL RESULTS 



To see the oscillatory behavior of the interparticle potential, we begin with the pure QED coupling case, i.e., we 
consider only Fig. 1(a) in F . In the following numerical calculations, only attractive interaction between two fermions 
with opposite sign is studied, so that Q1Q2 — — 47ra. 

Inserting the sum of Eqs.(|l]) and ([5]) into Eq.([T]), the pole of the integrand in the first quadrant can be identified 
numerically by replacing all q as q r + iqi in the denominator and then solving the complex equation. After that, 
one can calculate the undetermined parameters a and b through Eq. and consequently figure out the potential 
according to Eq. @. In Fig. 2 we demonstrated the real part q r (solid) and imaginary part % (dashed) of the pole as 
functions of temperature and the fermion mass. In FigJ^a), the mass has been fixed at 0.5MeV and in Figf^b) the 
temperature has been fixed at 0.2GeV. According to Eq.([2|), the fact q r and qi increase with temperature indicates 
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FIG. 2: Position of poles in terms of temperature and fermion mass in pure QED coupling, where the solid lines are for q r and 
the dashed lines are for §j. The curves in (a) are plotted with fixed fermion mass m — 0.5MeV and those in (b) are plotted 
with fixed temperature T = 0.2GeV. 
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FIG. 3: Pure QED interparticle potential oscillates in distance. The curves in (a) are for fixed mass m = 0.5MeV with the 
temperatures are 0.8GeV(dotted), 0.4GeV(dashed) and 0.2GeV(solid) respectively. The curves in (b) are for fixed temperature 
T = 0.2GeV with the effective fermion masses are 0.5MeV(solid), 0.8GeV(dashed) and 1.6GeV(dotted) respectively. 



the oscillation damps faster and more rapidly with increasing temperature. As to the fermion mass, there exist an 
opposite tendency, i.e., fast damping and rapid oscillations will be achieved by decreasing the mass. 

To make this analysis more visible, we choose three points on each curve in Fig[2]to plot the oscillatory potential 
in Fig|3J where FigJ3Ja) is fixed at m = 0.5MeV with the temperatures are 0.8GeV(dotted), 0.4GeV(dashed) and 
0.2GeV(solid) respectively. Fig|3jb) is fixed at T = 0.2GeV with the effective fermion masses are 0.5MeV(solid), 
0.8GeV(dashed) and 1.6GeV(dotted) respectively. Note that the solid curve in the two plots share the same set of 
parameters, and thus can be chosen as a reference. 

In Fig[31 one can see clearly the interparticle potential oscillates in distance, where the peaks and the oscillating 
damping tail vary with temperature and mass. The first peak, which is the most obvious one on the plot and can be 
regarded as the representation of all the peaks, decreases in amplitude and becomes broad in width with increasing 
temperature, which means the effective interacting distance becomes long though weak at relatively low temperature. 
Contrary to the temperature effect, the decreasing fermion mass enhances and narrows the peak as shown in Fig[3Jb). 

Now we are in the position to present the effect of bound states. To involve this effect, the contributions from the 
last two diagrams in Fig.l should be added to F, i.e., 

p _ pa , pa , pb+c , pb+c / *| -, \ 

vac 1 matt 1 vac 1 matt' V / 

Numerically figuring out the zero point of q 2 + F in Eq.([T]), one can obtain the oscillatory potential with respect 
to different group of parameters. In Fig|4l we exhibit q r (solid) and ^(dashed) in terms of T, a s , m and m s . The 
parameters we used in the calculation are presented below each plot. As a general impression, one can find out two 
tendencies: one is ascending with increasing parameters, like the first row in Figdl The other is descending with 
increasing parameters, like the second row. Based on the previous experience on pure QED coupling, one could expect 
the first tendency would enhance but narrow the peaks of interparticle potential and shift them to the left, which 
means strong but short distance interaction. While the second tendency would suppress but broaden the peaks and 
shift them to the right, which suggests the weak but long distance interaction. We show these tendencies more clearly 
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FIG. 4: Position of poles in terms of T, a s , m and m 3 , where the solid lines are for q r and the dashed lines for qi as in Fig[2] 
The parameters we used in the calculation are presented below each plot. 



in Fig[5]by directly exhibiting the potential in terms of these parameters where we set: a = y^y, m — lOMeV and 
T = 0.5GeV. 

One can observe that Figl5] is consistent with what we have argued by just analyzing the evolutional tendencies 
of the pole position with respect to various parameters. In each plot of Fig[5j the two curves denote for the case of 
rn s = 1.5GeV(solid) and m s = 0.3GcV(dashed) respectively, which shows the decreasing of m s or increasing of a s 
enhances but narrows the peaks on the potential and shift them to the left, indicating the strong but short distance 
interaction. Besides, it is worthwhile to compare the differences between the two curves in the same plot from the 
upper and lower ones which are mapped in different a s . The comparison shows the stronger the coupling a s is, the 
more the scalar mass can affect the shape of the potential, even the other parameters are the same. 

Finally let us compare the pure QED with the case involving scalar bound states. In Fig[Sl we demonstrated the 
potential at a s = (pure QED) and a s = 0.05 in upper plot with T = 0.3GeV and arrayed it with the lower one at 
T — 0.6GeV. The masses of the fermion and bound state are lOMeV and 0.4GeV respectively. There is no surprise 
that the pure QED potential is longer in the interactive range but weaker in amplitude than that involves the coupling 
with scalar bound state since it is the limit of the small coupling strength. 



IV. SUMMARY AND DISCUSSIONS 



We investigated the in-medium interparticle potential of gauge plasma with bound states by employing QED and 
sQED coupling, where complete one-loop boson self-energy has been involved which is equivalent to resum all one- 
loop diagrams. Different from Debye screening via HTL approximation, an oscillatory behavior of the interparticle 
potential with exponential damping amplitude has been demonstrated which is referred to the so-called Yukawa 
oscillation. The four parameters, a s , T, m, and m s in our model are free to change, since we do not intend to study 
the electromagnetic system but to employ it as a toy coupling form. In the numerical camputation we fixed a to the 
fine structure constant by considering the matter particle as point-like. While a s was treated as a effective parameter 
because of the inner structure of the scalar bound state. Evidently, a s — recovers the pure QED coupling. In the 
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FIG. 5: In-medium interparticle potential varies with the scalar mass at different a s . The solid line stands for the case of 
m s = 0.5GeV and the dotted line for the heavier mass case of m s = 1.5GeV. The difference between the upper and the lower 
frames lies in a s , where the former is for a s — 0.02 and the latter is for a s = 0.08. 



evaluation of the integral in Eq.([T]), we first found out the complex poles of the integrand, tracking their steps with 
parameters evolution. Then two opposite tendencies of the potential variation in terms of different parameters have 
been reported. The first tendency is to suppress the potential amplitude, slow down the oscillation "frequencies" and 
shift the peaks' centers to long distance. This tendency is caused by either decreasing a s and T or increasing m and 
m s . The second opposite tendency is to enhance the potential amplitude and frequencies, shifting the peaks center 
to short distance. As expected, this tendency is due to the increase of a s and T or decrease of m and m s . In brief, 
the potential between two fermions can be either long-weak or short-strong, according to the competition among the 
two groups of parameters. In some circumstances when the parameters are carefully "tuned", this potential is on the 
chance of being qualified strong and long enough to produce liquid state. 

We would like to point out here that the behaviors of the potential are governed by the static modes which 
are determined by the dispersion relation of plasma. The Debye screening, which describes the high temperature 
environment, can be obtained in the HTL approximation to the boson self-energy where the external momenta are 
soft. This implies only soft modes at q ~ gT get involved in the dielectric function when evaluating the interparticle 
potential. Therefore only those static modes in the soft region could be selected via dispersion relation and finally 
promise Debye screening. In other words, to reach Debye screening, one should take the high temperature limit 
by employing the HTL approximation from the beginning in the dielectric function to separate the soft modes, but 
not take the high temperature limit of q r and qi in Eq. f2| directly. As a matter of fact, the HTL static modes 
which are adapted to describe the extremely high temperature system may lose some significant details (hard mode 
contributions) in describing the medium at 1 ~ 3 times critical temperature. Instead, the improved calculation in 
this paper picks up the the long distance oscillation tail of the potential which might be helpful to understand the 
low viscous mechanism of QGP with scalar bound state. 

The interparticle potential is physically measurable at least in principle so we expect it is gauge independent. As 
to QED calculation, gauge problem does not emerge because the boson self-energy is gauge independent. While for 
QCD, the linear framework can only be persisted in the temporal axis gauge (TAG), but the one-loop calculations 
involving gluon self-coupling are various in different gauges even to the static modes. Anyway, we hope the discussion 
in TAG may give the qualitative features as well as correct directions of their evolution. More detailed investigations 
such as resummations are desirable in the future. 
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FIG. 6: Interparticle potential of QED and QED+sQED. The upper plot is at T = 0.3GeV and the lower plot is at T = 0.6GeV. 
The dotted curves in both plots are for the pure QED potential and the solid curves are for the potential involving scalar bound 
states where the coupling strength a s = 0.05. 
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